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CS| ■ Abstract 

rn ' 

^r^ ' We study a special kind of semiclassical limit of quantum dynamics on a circle and 

CN ' in a box (infinite potential well with hard walls) as the Planck constant tends to zero 

^" , and time tends to infinity. The results give detailed information about all stages of 

evolution of quantum wave packets: semiclassical motion, collapses, revivals, as well 
as intermediate stages. In particular, we rigorously justify the fact that the spatial 
distribution of a wave packet is most of the time close to uniform. This fact was 
previously known only from numerical calculations. 

We apply the obtained results to a problem of classical mechanics: deciding whether 

cd . recently suggested functional classical mechanics is preferable to traditional Newtonian 

one from the quantum-mechanical point of view. To do this, we study the semiclassical 

limit of the Husimi functions of quantum states. We show that functional mechanics 

remains valid at larger time scales than Newtonian one and, therefore, is preferable. 

Finally, we analyse the quantum dynamics in a box in case when the size of the box 
is known with a random error. We show that, in this case, the probability distribution 
of the position of a quantum particle is not almost periodic, but tends to a limit 
distribution as time indefinitely increases. 

1 Introduction 

In this paper we study some aspects of the dynamics of quantum wave packets in bounded 
domains. These themes are related to fundamental problems of theoretical and mathematical 
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physics. The dynamics of quantum systems in bounded domains has been studied for several 
decades [H 121 [3] and continues to attract attention [U O El [TJ [8]. Among the important 
results obtained are an analogue of the Poincare recurrence theorem for quantum systems 
with discrete energy spectrum and a detailed description of the structure of the revival 
phenomenon in such systems [21 El E], including the case of the infinite square potential 
well (box) m as well as the case of the finite square potential well j3] . Numerical simulation 
enables one to study in more detail the dynamical properties concerning collapses and revivals 
of quantum wave packets in bounded domains [6l [7] . Understanding quantum dynamics in 
bounded domains is important in condensed state physics and the physics of nanosystems 

There are some open problems in the quantum and classical dynamics of particles in 
bounded domains. In particular, it is pointed out in [6] that the issue of collapse of quantum 
wave packets has not yet been adequately studied. Another question: on which time inter- 
vals do the quantum and classical descriptions agree? The explicit form of the uncertainty 
relations for bounded domains is still an open question [TT]. The asymptotic properties of 
classical dynamics for collisionless continuous media in a box form the subject of papers 
by Poincare and Kozlov (see |l2l [IS HH HI]). Here we obtain analogues of the Kozlov's 
theorems on diffusion for quantum systems. 

We study the dynamics of quantum states on a circle and in a box using a special 
semiclassical limit as the Planck constant tends to zero and time tends to infinity (a similar 
procedure is used in the method of the stochastic limit P^). The results give detailed 
information about all stages of evolution of quantum wave packets: semiclassical motion, 
collapses, revivals, as well as intermediate stages. In particular, we rigorously justify the 
fact (previously known only from numerical calculations) that the spatial distribution of 
a wave packet is most of the time close to the uniform distribution (an analogue of the 
Kozlov's theorem on diffusion). This is done in Section |2] (the circle case, the main result 
being Theorem [1]) and Section [3] (the box case, the main result being Theorem [2|). We prove 
the theorems for coherent states on a circle and in a box and use the fact that an arbitrary 
wave function can be represented by an integral over coherent states. 

We then apply the obtained results to a problem of classical mechanics: deciding 
whether one should prefer recently suggested functional classical mechanics [T71 HB] (see 
also p^ [2D1 [211 [221 [231 [21]) to traditional Newtonian one. The basic concept of functional 
mechanics is not a material point or a trajectory but a probability density function in a phase 
space. Accordingly, the fundamental dynamical equations are not Newton (or, equivalently, 
Hamilton) equations but the Liouville equation (even if we consider just one particle, not an 
ensemble). The Newton (Hamilton) equations become approximate equations for the mean 
values of distributions of the positions and momenta. Corrections to solutions of the Newton 
equations have been calculated in some particular cases [T71 [HI [211 122] . 

Functional mechanics was suggested in an attempt to solve the irreversibility problem 
(or reversibility paradox), that is, to make the reversible microscopic dynamics compatible 
with the irreversible macroscopic dynamics (see [14j as well as ^5\ [25] )• This paradox is 
absent from functional mechanics since both the macro- and microscopic dynamical pictures 
become irreversible in some sense. 

A motivation of functional mechanics comes from the fact that arbitrary real numbers. 



being infinite decimals, are non-observable (and hence, so are individual trajectories). There- 
fore, it is more natural to consider bunches of trajectories (or the dynamics of the probability 
density) than individual trajectories of a material point. Each individual trajectory is a kind 
of "hidden variable" and has no direct physical meaning. 

A procedure of constructing the density function of a physical system from directly ob- 
servable quantities (results of measurements) is described in [23j. The dynamical interaction 
of the system and the measuring instrument is studied from the point of view of functional 
mechanics in [21] . 

In Section H] we try to approach the problem of choosing a preferable formulation of 
classical mechanics from a quantum-mechanical perspective. We ask whether Newtonian or 
functional classical dynamics remains consistent with quantum dynamics for longer time. To 
answer this, we study the semiclassical limit of the Husimi functions of quantum states of 
particles on a circle and in a box. For every h > a quantum density operator determines 
a classical density function on the phase space, and we pass to the limit as /i — )■ 0. Note 
that the evolution of the Wigner function and diffusion in collisionless media consisting of 
quantum particles on a non-compact space was considered in [26j . 

As a result, we obtain (Theorems [3] and H]) that both formulations of classical mechanics 
adequately describe the system when time is not arbitrarily large. But functional mechanics 
remains valid at larger time scale than traditional one. Hence, it is preferable in this aspect. 

Finally, in Section O we we analyse the quantum dynamics in a box in case when the 
size of the box is known with a random error (as we said before, we cannot know the exact 
size as an infinite decimal). We show that, in this case, the probability distribution of the 
position of a quantum particle is not almost periodic, but tends to a limit distribution as 
time indefinitely increases. 

2 Coherent states on a circle 

2.1 Definition of coherent states on a circle 

Consider a family of functions rjqp{x) G -^2(1^)5 (o'^p) G M^: 

1 \ {x — lY w{^ ~ 1) 



where a > 0, /i > 0. It satisfies a property known as the continuous resolution of unity 



"T // P[Vqp]dqdp= 1. 
^ri J Jr2 
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Here P[ip], ip G L2(M), stands for the one-dimensional operator acting on any vector 
(fi e L2(M) by the rule P[i']f = {ip, '^)4', where (■, ■) is the scalar product in L2(M). {P[4'] is 
a projector whenever ip is a unit vector.) The equality is understood in the weak sense: for 
all '?/', X e L2(M.) we have 

{tp, P[Vqp]x) dqdp = — — / / (^, Vqp){Vqp, X) dqdp = (^, x)- (2) 



2'Kh //e2'" "^^'^' ' ' 27ch 



In quantum mechanics, the functions T]qp (with fixed a) are called coherent states. The 
most general formal definition of coherent states was given by Klauder and Skagerstam [28] : 
a family of coherent states is defined as any family of vectors that continuously depend on 
their indices and form a resolution of unity. Another key feature of coherent states is that 
their properties are closest to those of classical particles among all pure quantum states (that 
is, all square- integrable functions). 

The following analogue of the family of coherent states for the spaces L2{—l,l) was 

introduced in |29j : 

+00 

^gp(^) = X] riqp{x-2nl), (3) 

n=— 00 

where q G [— /, /], p G R. 

Proposition 1. The functions ^ form a continuous resolution of unity in L2{—l,l): 



^ // P[vgp]dqdp= 1, 



27Th 

where Q = {{q,p)\q G [— /,/],p G M}. The equality is understood in the weak sense: for all 
'0, X G L2{—l,l) we have 

^ // ii^,P['Uqp]x)dqdp= ^ // {ip,Vgp){vgp,x)dqdp= (^,x)- (4) 

Here P[ip], ijj G L2(— /,/), stands for the one- dimensional operator acting on any vector 
ip G L2{—l,l) by the rule Plip]^ = {■ip,(p)ip, where (■,■) stands for the scalar product in 
L2HJ). 

Proof. We first establish the simple formula 

pi /'+00 

(^,^qp(2/-2m/)) = / ij{y)r]gp{y - 2ml)dy = / tp„iiy)Vqp{y)dy, 



where 

, , ]ip(x + 2ml), X e\-2ml -l,-2ml + 1], 

^ \0, x^ [2ml-l,2ml + l]. 

Using this formula and the property rigp{x — a) = rig+a,p{x) for every a G M, we get 

n n -|-00 n p 

{'4','^gp){vgp,x)dqdp= V] {ij,r]gp{y - 2{n + k)l)){r]qp{x - 2nl),x) dqdp 



X] I I {i^k, Tlq+2nLp)w{Tlq+2nl,p, Xo)]R dqdp 
k=—oo 

X] // ii^k,Vqph{Vqp,Xo)Rdqdp. 



k=—(:x) 



Here (■, ■)ir is the scalar product in L2(M). Then the desired equahty follows from 



+00 



{ilj,Vgp){vqp,x)dqdp= Y^ (^fc,Xo)R = (V^o,Xo)m = (^,x)- 

k=—oo 

The proposition is proved. D 

The functions Vqp may be called coherent states on a circle since, first, they satisfy the 
Klauder-Skagerstam general definition (continuous dependence on the indices and resolution 
of unity) , second, we shall see that the temporal evolution of these states on a circle tends 
to the dynamics of a classical particle on a circle in the semiclassical limit and, third, they 
converge to ordinary coherent states rjgp on a line as / — t- cxd. 

In Section [3] we show that the same properties hold for another family of functions in 

L2{-i,iy. 

+00 

These are coherent states in the infinite potential well (in a box). 

2.2 Spectral properties of coherent states on a circle 

One can express the functions Vgp in terms of the theta-function 

+00 
9{x,t) = 2. exp{— tttA;^ + 27ri/cx} 

fc=— 00 

(Re r > 0) as follows: 

{x — q)l pi P 



VZTra^ 
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J (x - g)2 ip{x - q) . 



The theta-function has the so-called modular property (Jacobi identity) [301 131] : 

9{-^-) = ^e^9{x,T). 
IT r 

Using this identity in ([5]), we obtain after some transformations that 

, . ijno^ f x-q pa^ 'Ka^\ f fap\'^\ 

M-) = V^^^-^p--^,^Jexp|-(-j j 

+00 ^ 2 1 

^ exp <^ -a^ (^-k - ^J + ijK^ - q) > 




(6) 



A;=— 00 



We have got an expansion of Vqp with respect to the orthonormal basis 



efc = ^e*t«^ A; = 0,±l,±2 
V2I 
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of the space L2{—l,l). The same result could be obtained by the direct calculation of the 
scalar product of Vgp with e^. For this purpose, note that 

I f + OO 

Vgp{x)e'^^'' dx = / r]qp{x)e'^^'' dx (7) 

■I J —00 



for integer k. Indeed, 

Vqp{x)e'^^'^ dx = ^ / r]gp{x - 2nl)e''T'''' dx 



+00 „_2nl+l p+00 

= ^ / r]gp{x)e''T'"' dx = / r]gp{x)e'T'''' dx. 

^^_^J-2nl-l J -00 

The integral on the right-hand side of ([7]) is a Gaussian integral, whose calculation also 
yields (jH]). Thus, the modular property of the theta- function turns out to be related to the 
Fourier transform of a periodic function defined by a sum of Gaussian functions. This is to 
be expected since the proof of the modular property involves the integration of the Gaussian 
function against trigonometric functions (see [30] and the more general case in [M])- In 
Appendix we give a proof of the modular property of the theta-function directly based on 
the Fourier series expansion of a periodic Gaussian function. 

The functions e^ are eigenfunctions of the self-adjoint operator 

TTC _ " " 



2m dx^ 
on L2{—l,l) with the domain 

D{H'^) = {^e AC\-l,l)\^{-l) = m, i^'i-l) = ^'(0), 

where ??i > is a constant. Here AC'^{—l,l) is the set of different iable functions whose 
derivatives lie in AC{—l,l), and AC{—l,l) is the set of absolutely continuous functions 
whose derivatives lie in L2(— /,/). The operator H'^ is a Hamiltonian (energy operator) 
for a free quantum particle of mass m on a circle (see [H]). Thus, the functions Vqp{x) 
can be expanded in a uniformly convergent series with respect to the eigenfunctions of the 
Hamiltonian operator of a free particle on a circle. 

Remark 1. By definition, a quantum particle on a circle cannot be free since there must be 
a potential constraining it to the circle. Only particles on the whole space (on the line in 
our one-dimensional setting) can be free. But we use the expression "free quantum particle 
on a circle" here to indicate the absence of potentials other than the constraining one. Also, 
as / — )■ 00, free motion on a circle becomes free motion on a line. 



2.3 General dynamical properties of the coherent states of a par- 
ticle on a circle 

The temporal evolution of the state Vqp is described by the formula 

'^qp,t = tVj Vqp, 

where 

U^ = eM-'{H^) (8) 

is the evolution operator for a free quantum particle on a circle. 

The function Vgp^t satisfies the Schrodinger equation with the periodic boundary condi- 
tions: 



ih- 



-^qp,t 10 U Uqp^i 



dt 2m dx"^ 

VgpA-^) = Vqp,t{l), <p,t(-0 = v'qpA^), 
Vqp,0{^) = Vqp{x), 

where x G [— /, /], t G M. Using the reflection method [32], we get 



(9) 



VqpA^) = X] '^IpA^ ~ '^^^) (^0) 



where 



''*'"' ry2.„ni + ny^ °'"'l- 4a.^(i + nO ^ H ^' '"' 

is a well-known function in L2(M) that describes the evolution of the initial wave packet rjqp. 
Here 7 = i^^- We see from (ITT!) that under free motion on a line, the centre of the wave 
packet moves along the classical trajectory q{t) = q + — while its dispersion grows, that 
is, the wave packet indefinitely spread with time. If Ag(0) = a is the initial mean square 
deviation of the coordinate, then its value at time t is equal to 



On the circle, only the terms with n = make an essentially non-zero contribution to 
the sum (fTO!) for small t (we assume that a <C /). Hence, we observe similar behaviour: the 
centre of the wave packet moves along the classical trajectory with period 

2/m 

P 

(the period of the motion of a classical particle of mass m and momentum p around the circle 
of circumference 21). The wave packet eventually collapses. The numerical experiments 
reported in [7] show that at time 

^ 2mla , , 

^"" - ^ '"' 

7 



one achieves an approximately uniform distribution of the position of a particle: 
\iJqp,t{^)\'^ ~ 1/2/. The value of Tcou is heuristically obtained as follows. The mean square 
deviation of the uniform spatial distribution on the circle [— /, /] is equal to -4=. Suppose that 

^Qii) = 7^ ) where Ag(t) is defined by ( IT2l) . Then t = Tcoii is a solution of this equation for 
t. Of course, this argument is non-rigorous since Ag(t) is the mean square deviation at time 
t for a particle on a line, not on a circle. But this conclusion is approximately confirmed by 
numerical calculations. We also obtain some rigorous asymptotic estimates corresponding to 
the flattening of the spatial density: see the next section. In particular, we will see that the 
distribution at time Tcoii (in the semiclassical approach) is non-uniform, the uniform density 
is achieved at slightly later times. Nevertheless, we will refer to Tcou as to a time scale of 
the wave packet collapse. 

The situation at large values of time is quite different from that of dynamics on a line. 
We express the solution of ([2]) as a series in eigenfunctions of H^: 



fc=— oo 



f .77 iht fTT ' 



The coefficients a^ ^p can be found from ([6]). We see that the dynamics is periodic: the wave 
packet is completely restored to its original form at the time 

_ 4m/^ 

that is, Vgp^Trev = '^qp- This phenomenon is referred to as the full revival of the wave packet. 
At time moments ^Trev with integer M and A^ one observes the so-called fractional revivals 
of the wave packet (|2l |3]): a copy of the original packet arises simultaneously at several 
places on the circle. (See [H 121 |3l IH El [7] for more details on the dynamics of quantum 
systems with discrete energy spectrum, including the particular case of quantum particles 
on a circle.) 

Thus, there are three time scales in the quantum dynamics of a particle on a circle [7]: 

1) Tci, the classical period of motion, 

2) Tcou, the characteristic time of collapse of the quantum wave packet, 

3) Trev, the full revival period of the quantum wave packet. 

In the next subsection we consider the semiclassical limit as/i— 7-0,0;— t-O, -— i-O (the 
parameter a occurs in the definition of Vgp] see dTj) and ([3])). The time scales have different 
asymptotic behaviour in this limit: 

Tcl = Ci, Tcoll = C2 — , Trev = ^— (15) 

n n 

where Ci, C2, C3 are constants. 



2.4 The semiclassical limit of the dynamics of coherent states on 
a circle 

In this subsection we consider limits in the space of distributions. Let K he a cyhnder (the 
set fl = [—1, /] X M 9 {g,p) with the points {—l,p) and {l,p) identified for all p G M). We are 
going to define the space of distributions on K. Introduce the set of test functions 

y{K) = {ff : M^ ^ M| 1) a{q + 2nl,p) = a{q,p), n = 0, ±1, ±2, . . . ; 

2)aeC^{R'y, 

3) ^^^ P \ si a so = '^' ^' *!' *2 = 0, 1, 2, . . .}. 

p-s>±oo dq^^dq-^2 



Its topology is defined by the seminorms 



P^(a) = max sup(l + ^2)^/2 

Si+S2<N m2 



d' 



1S1+S2 



a 



dq'^^dq 



S2 



N = 0,l,2,.... (17) 



This is the space of functions on K that decay rapidly with respect to p. Let J^'{K) be the 
space of distributions, that is, continuous linear functionals on ^{K). 
Let a G S^{K). We define the following distributions: 

iS{q-qo,p-po),cr) = cr{qo,po), 

/+00 
fiq)a{q,po), 
-00 

{cS{p-po),(T) = c a{q,po)dq, 



where {qo,po) G K , /(g) is an integrable function on a line, and c G M. Consider the function 
where D G (0, 00), and define it for D = and D = 00 by putting fo{q) = lim foiQ) = ^(?) 

D— >0 

and ^Poo{(l) = hm ^d{(1) = h (the limits are taken in the space y'{K) of distributions). 

The space S!{K) of test functions is defined by the same formula as y{K) but with 
the third condition of rapid decay with respect to p in f lTB]) replaced by the condition of 
being compactly supported with respect to p. We similarly introduce the space ^'{K) of 
distributions. 

Theorem 1. We have the following limit formula in S^'{K) (where {q,p) are fixed and 
{q',p') are variables of integration with test functions cr{q',p') G J>^{K)): 

1 1 ^^-^^ 2kl n 

lim{^IKp,^g'p',i)P-^ Yl VdW -q-j^-a+—{t-cTrev)W -p)} = 0. (19) 

fc=0 



The limit is performed as follows: /?. — )■ 0, a — )■ 0, - — )■ 0, t = t{h), -(t — cTrev) -^ 2mD, 
h{t — ^ ) — )■ 0, where c G M, D & [0, oo], and the numbers N' and a depend on c. If c is 
rational (and hence, can be written as a reduced fraction c = ^), then N' = N for odd N 
and N' = Y for even N. Further, a = ^ for N = 2 (mod 4) and a = otherwise. If c is 
irrational, then N' = 1, a = 0. The parameter a occurs in the definition ofvgp (see formulae 
(Cy and ^), Trev = ^^ (see p4^)- The convergence in 17^) is uniform with respect to 

{q,p)en. 

Let us give some comments on Theorem [H The quantity 2^|('^gp5'^g'p',t)P is the prob- 
abihty density for a quantum particle to be in the state Vqp at time t under the condition 
that it was in the state Vq/p/ at time (up to the norms H^gplP and H'Uij'p'll^, which tend to 1 
in our hmit; see Proposition |2] below). 

The limit /i — ;■ corresponds to the semiclassical approximation. The limits a — )■ 0, 
- — ;• correspond to convergence of the mean square deviations of the position and momen- 
tum of the quantum wave packet Vqp to zero. We actually have Ag ~ a and Ap ~ ^, where 
Ag and Ap are the mean square deviations of the position and momentum (respectively) for 
Vqp for all (g,p) G fi (see [H]). Here the notation f ^ g means that lim - = 1. Thus, in the 
semiclassical limit under consideration, a quantum particle in the state Vqp has a well-defined 
position (equal to q) and momentum (equal to p) just as classical particles do. Therefore, 
we may say that, in the semiclassical limit, ^^K'^gp, ^(j'p',t)P is the probability density for a 
quantum particle on a circle to be at the phase point (g,p) at time t under the condition 
that it was at the phase point {q',p') at time 0. 

In Theorem [1] we consider various rates of convergence of t to infinity with respect to h 
and a, that is, the various time scales l)-3) listed in the end of Subsection 12. 3[ 

We see from (fTS!) that the case when c = and D = (— — )■ 0) corresponds to the 
classical time scale T^;: time is either fixed or increases slower than the decrease of the 
collapse velocity (proportional to -) of the packet. Then formula 019p takes the form 

1 T) 

lim[— — I (t;gp, t;g'p',t)P - ^{(l' -q + —t,p'- p)] = 0. 

Here 5{q' — q + —,p' — p) is the probability density for a classical particle on a circle to be 
at the phase point (g,p) at time t under the condition that it was at the phase point {q',p') 
at time 0. Thus, in the semiclassical limit at time scale Td, we have classical dynamics: the 
quantum probability density of transition to the phase point (g, p) for a particle that was at 
the phase point {q',p') at time is equal to the corresponding classical probability density. 
The case when c = and D E (0, oo) (— -> 2mD G (0, oo)) corresponds to the time 
scale Tcou- Formula flT^ takes the form 

1 p 

lim[^|(t^gp,t^g'p',t)P - ^niq' ~ 1 + —^) ^{p' - P)] = 0- 

We see that the quantum probability density of transition from one point to another is 
already different from the classical one: there is a spatial spread of probability distribution. 
In particular, the case t = Tcoii corresponds to D = -j= and, since D is related to the mean 
square deviation of the "spreading function" ^j^, this deviation becomes approximately equal 

10 



to the mean square deviation of the uniform distribution on [—/,/], which agrees with the 
numerical results in [6l [^ . However, the semiclassical limit fails to provide the exact uniform 
distribution, in contrast to the following case. 

The case when c = and D = oo (— — )■ oo) corresponds to an intermediate time scale 
between Tcou and T^ev Formula (HM takes the form 



1 



l^TTT^IKp'^a'p',*)!^ = 7T7^b'~P)' (20) 



so that Theorem [T] corresponds to a complete flattening of the spatial probability density 
in this case. Therefore, we also associate this case with time scale Tcou (corresponding to 
collapse of the localized wave packet). Note that here we get a mathematical justification 
of the (asymptotic) flattening of the spatial probability density for a quantum particle in 
a finite volume, which was previously known only from numerical calculations [6l [7]. This 
result may be regarded as a quantum analogue of the Kozlov's theorems on diffusion for 
classical systems [131 [IH US] • 

The case c 7^ corresponds to the time scale T^ev If c is irrational, then (just as in the 
previous case) formula (TT^ reduces to (12U1) . that is, one observes a complete flattening of the 
spatial probability density. The case of rational c corresponds to a revival (fractional or full) 



of the wave packet. We discuss this case in more detail. First suppose that t — -rrTrev — ^ 



N- 



whence D = (in the simplest case t = ^Trev)- Then formula flT9l) takes the form 



N'~l 



lim 



\{vqp,VgY,t)\'^ = — Y^ S{q' -q- — -a,p' -p), 



A:=0 

where N' and a are defined in Theorem ^ N' = 1 (that is, c is integer or half-integer) 
corresponds to a full revival of the packet, A^' > 1 corresponds to a fractional revival of the 
packet. This agrees with the results in [2l [3l Hj. In the case D G (0, 00) we obtain a sum of 
A^' "spread" wave packets (see (IT^ . the formula is not simphfied in this case): the revived 
wave packets begin to spread. In the case D = 00 we again get a complete flattening of the 
spatial distribution density (120]) . 

Since every irrational number can be approximated by rationals within any accuracy, we 
can (loosely) say that the case of irrational c in Theorem [1] is a limiting case of rational 
c as A^ — )■ 00: if we approximate an irrational number by a sequence of rationals, then 
their denominators increase, the distance between neighbouring terms in the sum of delta- 
functions in flT^ tends to zero, and the sum of the delta-functions tends to the uniform 
distribution (in the weak sense). In other words, the cases of an irrational c and a very close 
rational c' are almost indistinguishable. 

Note that the distribution of the momentum is preserved in all these limiting cases. 

Thus, we have traced the whole evolution of a quantum wave packet on a circle. A well- 
localized initial wave packet eventually collapses until there is a complete flattening of the 
spatial density. At certain moments we see that copies of the initial packet simultaneously 
arise at several points of the circle and then again eventually collapse until there is a complete 
flattening of the density. Since there are "more" irrationals than rationals, we can say that 
the particle most often stays in states whose spatial distribution is close to uniform. 

11 



Thus, Theorem [T] completely describes the free quantum dynamics of a particle on a 
circle at all time scales in semi classical limit. All stages are parametrized by the two real 
parameters c and D. 

Here is a simplified version of Theorem [H which deals only with principal time scales 
(the classical motion, complete flattening, and exact revivals) without intermediate ones. 

Corollary 1. We have the following limit formulae in S^'{K): 

1) 

X T) 

iim[^IKp,t;gv,t)|^ - Kq -q + —'t,p'- p)] = 



as ^, a, - — 7- 0, t = const; 

2) 



lim^lKp,t;gV,t)|^ = Yj^^ijp' -V) 



as h, a,- ^ 0, t ^ oo, — — ;■ cx), 

as well as h,a,- -^ 0, t = cTrev —^ oo, where c is irrational and Trev = ^^^ 
3) 



^'' ir77 5Z^(^'-^-^-«'P'-P)]=0 



k=0 



as h,a,^ ^ 0, t = cTrev -^ oo, where c = j^ is rational and the numbers N' and a depend 
on N. 

These limits are uniform with respect to {q,p) G Q. 

Proof of TheoremUl Let us express the scalar product of Vgp and Vg'p>^t in L2{—1, 1) in terms 
of the scalar products of rjgp and rigi^2ki,p',t, k = 0, ±1, ±2, . . ., in L2(M): 

{vqp,Vq'p\t) = ^ / ri^{x-2kl)r]qipi^t{x-2{k + n)l)dx 



= X^ / V^{.x)riq>j^2kl,p',t{.^)dx = y^ iVqp,Vq'+2kl,p',t) 

fc=— oo ~°° fc=— oo 

Substituting in (J2T1) the expression for the scalar product 
{Vqp,Vq'p',t) = X —-—exp > 



(21) 



2 + i7 "^ 1 4a2(2 + z7) 2%^ 

iip' + p)iq' - q) it{p' + pY 



2h 8mh 



(22) 



where 7 = ^-^ as above, we get 



i{p' +p){q' -q + 2kl) it{p' + pf 
2h 8mh 
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['^qp^ Vq'p\t 



,2_ 2 g ^^pj W-g + 2kl+^J^^^ 



v/4T^ 



fc,n=— oo 



4a2(2 + z7) 



{q'-q + 2nl + ^^^f a\v' - vf ^{p' +p){k- n)l 



4a2(2 - 27) ;i2 ;^ 

Put k = n + r and replace the sum over k and n by a sum over r and n: 



[Vqp, Vq'pi^t 



+00 



exp 



{q' - q + 2nl + rl + 



{p'+p)t \2 
2m > 



V ' ' r,n=— 00 k. \ ' y 

r2/2 iH7(g' - g + 2n/ + H + ^^^) a^p' - p)^ i{p' + p)rl 



a2(4 + 72) 



^2(4 + 72^ 



^2 



;i 



Let cr{q',p') be an arbitrary test function from ^{K). We expand it in a Fourier series 
with respect to q' and represent it by a Fourier integral with respect to p': 



+°o /•+00 



1 J^^ /•+00 

(^{q', p) = —r^ Y] / «j (^) exp { i-jq' + iz/p' I di/. 



j=-oo 



(23) 



We calculate the following integral: 



{vqp, Vq'p'^t)\'^ exp I i-jq + zi/p 



' ^ dq'dp 



n 



[vqp,v^,_i^,,,)rexpyjj [<l-^)+^P ( -- 2;^ 



pt 



njt 



dq'dp' 



+00 



exp 



(g' — g + r/)' 



r^/^ irljlq' — q + rl) 



a^ip' — p)"^ i{p' + p)rl ,TC 



^2 
+00 
27r/l >^ exp 



n 



+¥l^'-|^)+^^r-i^)^^^'^^' 



«2(4 + 72) 



r/7 



+ 



TTJ 



H 2 



a2(4 + 72) / 



4a2 



, TTJht , 

2lm 



1 2 



2iprl 



2/2 



r^/ 



+ zyj (g-r/-— ) +zpz/J>. (24) 
a^(4 + 7^) f V "^ 



We now realize all the passages to the limit. First consider the case when c = 0, 
D G (0, 00). Only the term with r = remains non-zero in fl24l) . We have 



lim 



^ / / I (^gp> ^g'p'.t) P exp j^yjg' + * V| f^^'^^P' 



-exp^-l(^D 



^yJ [Q 1] + ipv 
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This means that 
1 



hm 



27r/i././f^ 



fgp, fg'p',i)l (y{q\p')dq'dp 

+ 00 „_|_oo 



L= ^ |_^ a,(.) exp -i (^d) + .^, (^ - ^) + ^P- d. 



2v^ 



0. 



Clearly, the integrals and the series converge uniformly in (g,p) G ^2. Again expressing aj{v) 
in terms of cr(g',p') by the formula 



aj{v} 



-1= jj (T{q',p')expi^-ijjq' -iup'jdq'dp' 



and using the modular property of the theta-function, we obtain (TT^ . A similar and simpler 
argument proves flTOj) for D = and D = 00 (here c = as above) as well as for irrational c. 
Consider the case of an arbitrary rational c = ^. We first assume for simplicity that 
t — ^Trev -^ (whence D = 0). Then the terms of 1^^ with non-zero limits are only those 
with r = —2cj (this can be seen from the two first terms in the exponent). Accordingly, the 
terms of the sum ( l23l) have zero limits unless j is such that 2cj is an integer. Namely, the 
terms with non-zero limits are those with j = N'J, J = 0, ±1, ±2, . . ., where N' = N for 
odd A^ and N' = ^ for even A^. We look at the term —injr = 27iicj'^ in the exponent of 
(12^1) in more detail. If A^ is odd, then exp{2nicf) = exp{27riMNJ'^) = 1 for all J. If A^ is 
even, then exp(27ricj^) = exp(i7rMA^'J^) = (— 1)"^ = exp(27rA^'J) (if A^ is even, M must be 
odd since M and A^ are coprime). Thus, for odd A^ we get 

+'^ p+00 



^^T^^^ \i^qp,Vq'p',t)\'^'^{q',p')dq'dp' = —== Y^ / aNj{j^)exp^ijNJq + ipuy 

(25) 



j=— 00 ' 



This is equivalent to saying that 



N-l 



lim^-^ // \(^<iP^^<i'p',t)\'^^iq'^P')dq'dp' = —J2^i(l+ ^^P) 



2nh ki V-. / . . N ^ '' N 



This proves formula flT9|) for odd A^. If A^ is even, we get 



+00 
i 

2v^ , 
This is equivalent to the equation 

Af'-l 



^ +00 ^+00 ^ 

-7= Yl / aN'Ai^)exp\i-N'J{q + l) + ipu\ . (26) 



J J it' J — n 



fc=0 
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If A^' is even (so that A^ is divisible by 4), then / + ^ = for some k. If A^' is odd (so that 
N = 2 (mod 4)), then / + ^ = |^ for some k. Therefore, we can write 

N' — l 

where a = if A^ is divisible by 4 and a = ||- if A^ = 2 (mod 4). This proves formula (I19p 
for even A^. 

The same argument works when the condition t — ^Trev -^ does not hold, but still 
we have D = 0. The only difference is that the first argument of a acquires an additional 
summand ^{t— ^T^ev) in the final expression (because the sum ^^ + ^^^ in the exponent 
of (1211) does not tend to zero in this case). 

We similarly treat the case when D G (0, oo). Again, the terms of (I24p (resp. of the sum 
( 123|) ) have zero limits unless r = —2cj (resp. j is such that 2cj is an integer). However, just 
as in the case when c = 0, the integrands of the right-hand sides of (p5|) and ( |26|1 acquire 
a factor e~2^~r^) ^ which results in the replacement of the delta-functions by the functions 
ifD- In the case when D = oo, all terms of the sum over k tend to ^5(p' — p), which yields 
formula ( TT9|) for the last limiting case. 

The theorem is proved. D 

We now prove another proposition to be used in what follows. 
Proposition 2. The norm of Vgp tends to unity uniformly on Vl as h ^ 0, a —^ 0, - — > 0. 
Proof. By formula f l2T]) we have 

fe=— oo ^ ^ 

n 

Here are some remarks on Theorem [H 

Remark 2. We refer to the limit performed as to the special semiclassical limit. The usual 
well-known semiclassical limit (see, for example, |33]) is the limit ^ — )■ with constant t. 
In our case, this corresponds to the first time scale Td (classical motion). We consider the 
simultaneous limits h -^ and t — > oo, where t and h are related to each other and to 
a — )■ in certain different ways. This allowed us to investigate analytically not only the first 
(classical) time scale, but the other two time scales and the intermediate time scales as well. 

Remark 3. We have proved the theorem on the semiclassical limit of dynamics on a circle 
only for the quantum states of special form. But the result obtained can be extended to 
more general states since, by formula (jl]), every function ip G L2{—l,l) can be represented 
by an integral over coherent states: 

^ = ^ / / M(i^p)^qpdqdp, 
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where U{q,p) = Kp,^). 

Nevertheless, it is worthwhile to give a formulation and a proof of Theorem [1] which would 
not distinguish any particular form of wave packets at all|j This is a subject for further work. 

Remark 4. Consider the dynamics of the mean position and mean momentum of the state 
Vqp^f The momentum is easily seen to be preserved under our semi classical limit /i, a, ^ — ?■ 
at all time scales, that is, for any behaviour of the time variable t: 

p = -ih / Vgp^t{x)v'gpt{x) dx -)■ p. 

We now look at the dynamics of the mean position. Since the position operator is not 
well-defined for a particle on a circle, one usually considers the mean complex exponent of 
the position, which is uniquely determined: 



(^'^'), 



e'i''\vqpAx)?dx. 



Consider the time scale T^, that is, -^ — )■ 0. We assume for simplicity that t = const. Since 
l^qp^tix)]"^ tends to the periodic delta-function J2n ^(^ ~ 1 ~ m^ ~ ^^0 ^^ this time scale, we 
have 

Thus, at time scale T^, the centre of the wave packet moves along the classical trajectory. 
At other times scales, the notions of "centre of the packet" and "mean position" have no 
physical meaning since the packet either collapses and has no definite centre or has several 
centres (fractional revivals). In the last case the usual expectation of the (complex exponent 
of the) position may lie between these centres. This also lacks physical meaning because the 
particle cannot be observed near this expected value. Only a full revival makes the notion 
of position meaningful again. 

Remark 5. We use the limit /i, a, - — )■ as a mathematical tool. In particular, it allows us 
to formulate the result about the flattening of the spatial density. As already mentioned, 
the fact that the spatial density of a particle in a finite volume is most of the time close 
to uniform, has not yet been proved in the literature nor even stated in a mathematically 
rigorous manner: the uniform distribution is never achieved exactly for a localized initial 
packet. What then is the meaning of the "closeness" of our distribution to the uniform one? 
We have proved that the spatial density is exactly uniform in our limiting case at certain 
time scales (scale Tcou and most of scale Trev)- This is a mathematical expression of the fact 
that the spatial density is close to uniform. 

These limits have no direct meaning from the physical point of view: Planck constant h is 
a physical constant and cannot tend to zero (see [31] for a discussion of the semiclassical limit 
with the constant h for the baker's map). The parameter a is also fixed for a fixed coherent 
state. From a physical point of view, it would be more correct to consider limiting cases 
for dimensionless quantities [35]. Let us restate our results in terms of relations between 



^The authors are grateful to J. R. Klauder for this important remark. 

16 



the time scales Td, Tcou, T^ev and time parameter t. Using formula ( !T5|) . we can rewrite the 



limits /i — ;• 0, a — )■ 0, - — )• in the form 



-^ coll ^ rev ^ rev 

-— > OO, > CXD, — > OO, 

J-rl -L ml! -L rl 



t d -f coll 

or 

-tret) ^^ -'-coll ^ -^ ch 

that is, the time scales (the classical period of motion, characteristic time of collapse, 

and full revival time) are distant from each other on the time line. Moreover, the limits 

^(t — cTrev) — ^ 2mD and ht(t — cTrev) — ^ occurring in Theorem [1] can be rewritten in the 

form 

f — cT f 

V (^-i ret) i-> 7-^ 

> 2mD, > c. 

-t coll -t rev 

Thus, we see that the limits under consideration have a clear physical meaning. 

3 Coherent states in the infinite well 

3.1 A map of the dynamics of a particle moving on a circle to the 
dynamics of a particle moving in a box (classical mechanics) 

Consider a particle moving freely in a one-dimensional infinite potential well [— /, /] with rigid 
(elastic) walls. In classical mechanics, this case can be reduced to dynamics on a circle. To 
do this, one uses a two-sheeted covering of an interval by a circle (see, for example, [131 113] )■ 
If we are given an interval [— /, /] with phase variables {q,p) and a circle [—21, 21] (the points 
—21 and 2/ being identified) with phase variables {q',p'), then the two-sheeted covering of 
the interval by the circle is determined by the formula 

p'=\p\. (28) 

Here we first shift the interval [—1, 1] to [—21, 0], then reflect it in the point and obtain the 
interval [—2/, 2/], and then glue the points —21 and 21. Thus, instead of oscillations on the 
interval, we get a rotation in one direction on a circle. We have reduced the dynamics of 
a particle in a box to that of a particle on a circle. However, we note that the sign of the 
momentum is not uniquely determined at g = it/ since it jumps at these points. This is a 
consequence of the approximate nature of the model. 

To define the spaces of test functions and distributions on the strip Q, we first introduce 
the space of rapidly decaying functions on fl: 

y{n) = {(t:M2^M| l)(T[(-l)"(g + 2n/), (-!)>] = a{q,p),n = 0,±1, ±2, .. .; 

2)ae C°°(M2) 




d- 



■S1+S2 



(29) 



a 



3) ^^^ P \ s.fi so = 0' ^' *i' ^2 = 0, 1,2,.. .}. 
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The topology on this space is given by seminorms (TT7|) . The space S^'{Q) of distributions is 
the space of continuous hnear functionals on ^{fl). 
Let a G y{Q). We define the following distributions: 

{6{q-qo,p-po),a) = a{qo,Po), 

r+oo 

{f{Q)S{p-po),(r) = fiq)(r{q,Po), 



J — oo 

{cS{p-Po),(t) = c a{q,po)dq, 

where {qo,Po) ^ I^^? Hq) is an integrable function on a line, and c G M. 
We again consider the function 

where D G (0, oo), and define the distribution ^PD{q)S{p — Po), Po & ^, at D = and D = oo 
using the corresponding limits in the new space of distributions: 

M(l)S{p-Po) = lini¥'D(g) = S{q,p-po), 

fooiq) = lim fniq) = -rjKP - Po) + tAp + Po) 

(the limits are taken in y'{Vi)). 

Functions in y{VL) have the following important property: 

a{±l,p) = a{±l,-p). (30) 

Let K21 be the cylinder K (defined in Subsection 12. 4p with / replaced by 21. In other 
words, K21 is the set ^^2/ = [—2/, 21] x M with the points {21, p) and (— 2/,p) identified for all 
p G M. The space ^(i^2/) is just the space y{K) (defined in ( !T6|) ) with / replaced by 21. 
We easily see that S^iVt) is a subset of y{K2i). Namely, y{Q) = Ty{K2i), where the map 
T is defined by the formula 

a{q,p) = T[a]{q,p) = p{q ~l,p)+ p{l - q, -p). (31) 

We define 

(' '\ T--ir u ' '\ /|^(^ + ^''P'). g' < 0, 

y2^i}-q,-p), g' > 0. 

Here p G ^(7^2/), o" G ^{Q). The map T^^ is not an inverse: we have TT~^a = a 
for all functions o" on a circle but T~^Tp = p only for even functions, that is, those with 

p{(i',p') = p{-q',-p')- 

We also introduce the space of test functions ^{Q) by the same formula fl29l) as for ^(Q) 
but with the third condition of rapid decay with respect to p replaced by the condition of be- 
ing compactly supported with respect to p. We similarly introduce the space of distributions 
^'(f^). We also have ^(fi) = T^{K2i). 
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3.2 A map of the dynamics of a particle moving on a circle to the 
dynamics of a particle moving in a box (quantum mechanics) 

We note that finding the position (on a circle) of an image under the map fl28|) requires the 
simultaneous knowledge of the position q and momentum p (or rather, the direction of the 
momentum) in a box, which are, thus, presupposed to be simultaneously well defined. But 
this is known not to be the case in quantum theory by the uncertainty relations. There- 
fore, one cannot use the map fl2Sl) directly to construct the corresponding map in quantum 
mechanics. 

We define another map 

e : L2(-2/, 21) ^ L2(-/, /), V^(x) ^ [e^](i/) = ^[^(y - /) - ^{i - y)\ 



and a map 



^(^(x + O, x<0. 



<$>-': u{-i^)^u{-ii,ii), ^(2/)^[e-V](^) = <_'^^^^;_;^ ^-q' 

The map B is an analogue of T and G~^ is an analogue of T~^ (see (13T]) and (!32|) ). 
Clearly, we have 6B~^(y9 = y9 for all functions </) G L2(— /,/) but Q^^O?/' = i\} only for odd 
functions tp G L2{— 21,21). Restricting G to the set of all odd functions in L2{— 21,21), we 
get a one-to-one correspondence between odd functions in L2{— 21,21) and all functions in 
L2{—l,l), which is just the usual odd extension of functions to the interval of twice the 
original length. However, we define B for all functions from L2{— 21,21). 

Proposition 3. 1) The map G~^ preserves scalar products: 

{<c)~^ip,<c)'^>c)2i = {^,>t:)i 

for all ip,>c G L2(— /,/). Here the subscripts 21 and I indicate the scalar products in the 
Hilhert spaces L2{— 21,21) and L2{—l,l), correspondingly. 
2) The map 6 preserves scalar products: 

(ev',ex)/ = (^,x)2z, 

if at least one of the functions ip G L2{—21, 21) or x ^ L2{—21, 21) is odd. 
Proof. 1) We have 

[Q-M{x)[Q-^K\{x)dx 

1 r - 1 /"''- 

= - / {p{x + l)x[x + /) dx H — / (p{l — x)h{1 — x) dx 
2 J -21 2 Jo 

Tp{y)>i{y)dy= {^,k)i; 
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2) For definiteness, assume that ip is odd. Then 



= \j \^{y - - V^(^ - y)My - - x{i - y)]dy 

i'iy - l)x{y - Ody + / i^il - y)x{l - y)dy 
' -I J-i 

ipix)xix)dx = {ip,x)2i- 

-21 



n 



3.3 Coherent states in a box 

Consider the coherent states Vqp (defined in (3)) on the space L2{— 21,21): 

+00 
^gp(^) = X^ r]qp{x-Anl), 

n=—oo 

where {q,p) G ^22^. Unless otherwise stated, we always assume in this subsection that Vqp 
are the functions from L2{—21, 21) defined by ([3]) with / replaced by 21. 
Then we easily see that 



-^u;/_g,_p(y), g>0 



PVgp\iy) = \ V2 . /,.N „. . n (33) 



(we have used the property rjqp{x) = ri^q^^p{—x)). Here 



+00 



^qp 



(y) = E I'^^P^y - 4^0 - Vqpi2l - y + 4n/)] = J2 (-l)%p[(-l)"(y - 2nl)], (34) 



for {q,p) G Q. One can similarly prove that 



,-1 n. ^ V2, 



[6 Uqp]{x) = —[Vq-l,p{x) - Vl-g,^p{x)] (35) 

(that is, <d~^Qvqp 7^ Vqp since the function Vqp is not odd). Note that ^±^,0 = 0. 

We now prove an analogue of Proposition [1] (which actually follows from that proposition 
and properties of the map 9). 

Proposition 4. The family of functions [34\ ) forms a continuous resolution of unity in 
L2i-lJ): 



— // P[ujqp]dqdp = 1. 



27ih^ _^ 
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Equality here is understood in the weak sense: for all </), x G L2{—1, 1) we have 

^ / / (V'' P[^qp]>^) dqdp = ^ / / {(fi, ujqp){ujgp, k) dqdp = {if, x). (36) 

Proof. Using the properties of G, O^^ and Proposition [H we have 
{(p, x) = (0"V, Q'^>i) = ^ / / (0"V, ^gp)Kp, e"V) dqdp = 



2Txh 



{ip,Qvqp){eVgp, x) dqdp 



AttH 



(V3, Ug+i^p){ujg+i^p, x) dqdp + 



A-nh 



{(p,UJl^q^^p){uJl^ 



q,~p^ 



X 



) dqdp 



2'iih.,.,Q 



{ip,ujqp){ujgp,x)dqdp, 

as required. Here we have written Q^i = [~2/, 0] x R and Q21 — [0) 2/] x M. D 

Proposition|l]along with some properties concerning semiclassical dynamics (to be proved 
below; see Theorem [2]) enables us to call the functions Ugp, {q,p) G fi, coherent states in the 
infinite potential well (in a box). 

One can obtain an analogue of formula (^ for the functions Ugp-. 

{x-q)l 2pl 4/2 



UJgp{x) 



</2 



vra^ 



^ 



e 



■.e 



ma' 



nh^ vra^ 



exp 



vra^ 



{x-2l + q)l 2pl 4/ 



2 1 



ma^ 



Tch ' vra^ 



exp 



(x — qY ip{x — q) 
Aa^ h 

{x-2l + g)2 ip{x -2l + q) 



Aa^ 



h 



Using the modular property of the theta-function, we get the following analogue of for- 
mula (EI): 



UJ. 



qp 



7ia^ 
32/^ 



X — q pa^ 'Ko? 



Al 



2ilfi AP 



X — 21 + q pa^ Tta^ 




+00 

E 

fe=— 00 
— exp 



exp^-a^(-fc-^ 



vr 



Al 2ilh' AP 

p\2 i7[k{q — I) 



-{^? 



21 



/ TT p\ 2 iTTk{q — I) 



" l2/' + 



h 



+ 



21 



sm|^(x-/)). (37) 



This yields an expansion of ujgp with respect to the orthogonal basis 



^'= = 7i'^"(ir(^-^) 



1,2 



As in the circle case, the same result can be obtained by a direct calculation of the scalar 
product of ujgp with the elements of this basis if we use the formula 

/ ujgp{x) sin ( — (x - /) j rfx = / rigp{x) sin {—{x-l)\dx 
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for integer k. This formula is obtained from the expression 

Vqp{x)sm(—kx]dx= / ■r]qp{x) sin i—kx) dx 
(which is proved in a similar way to ([7])) using properties of 0: 

/ r]qp{x) smi—{x-l)jdx= / Vg-iA^) sin f -^x j dx 
= / Vq-i^x) sin i^kx) dx = fsin (-^kx) ,Vq-iA 
sin y-^^kxj , Qvg^i^pJ = (^sin y^i^ix - I)) , uj, 



ujgp{x)sm ( ^(^ ^0 ) dx. 



e 

x;gp(x) sm ( 
The functions /^ are eigenf unctions of the operator 



9P, 



2m dx^ 

with the domain 

D{H^) = {^e AC\~l, /)| V(-0 = ^-(0 = 0}. 

The operator H^ is a Hamiltonian (energy operator) for a free (see Remark [1] for the use 
of the notion "free") quantum particle in the infinite square potential well (see [H]). Thus, 
the functions u)qp{x) can be expanded in uniformly convergent series with respect to the 
eigenfunctions of the Hamiltonian for a particle in a box. The evolution of the state Uqp over 
time is given by the formula 

^qp,t = U^ Uqp, 

where 

U^ = eM~H') (38) 

is the evolution operator for a free quantum particle in a box. The function ujqp^t satisfies 
the Schrodinger equation with the boundary conditions corresponding to the infinite square 
well: 

■ ^dujqp^t _ ha UJqp^t 

dt 2m dx'^ 

(^qpA-^) = ^qp,S) =0) 
^qpfli^) = Wgp(x), 

where x G [— /, /], t G M. Using the reflection method, we arrive at 



U^qpA^)= E (-l)%p,*[(-l)"(^-2ri/)]. 



22 



where rjqp^t is defined as above (see (ITT!) ). 

Note that formulae (133|) and ( 135|) remain valid if we replace t;^^ and Uqp by t;gp,t and Uqp^t 
(respectively) for an arbitrary t. 

Expanding in eigenf unctions of H^, we get 



(TT \ I 1 Tit I TT K" 

yA;xjexp<^-— ( — 



2 



where the coefficients bk,qp can be found from formula (1371) . 

The dynamics of coherent states Uqp^t in a box is similar to that of coherent states on a 
circle. For small values of time, a well-localized initial wave packet remains well-localized 
and its centre moves along the classical trajectory with period 

Aim 

J-d — • 

P 

The wave packet eventually collapses, and at the moment 

2mla 



Z 



coll 



V3h 

we observe an approximately uniform spatial distribution. However, at the moment 



T 



irh 



we observe a full revival of the wave packet: Uqp^Trev = ^gp- Fractional revivals are observed 



at the moments ^Trev for all integer M and A^. The full revival time now increases by a 
factor of 4 compared to that for motion on a circle (compare the last formula for Trev with 
formula ( !T4|) ) because the energy spectrum is now twice as dense as that for a particle on 
the circle [— /, /]. 

Thus, as in the case of a particle on a circle, we have three time scales: 

1) Tci, the classical period of motion; 

2) Tcoii, the characteristic time of collapse of the quantum wave packet; 

3) Trev, the period of full revival of the quantum wave packet. 

The asymptotic behaviour of these time scales as /i, a, - — )■ is again described by formula 
([15]). 



3.4 The semiclassical limit of the dynamics of coherent states of 
a particle in a box 

We now prove an analogue of Theorem [1] for the box. 
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Theorem 2. We have the following limit formula in y'{VL) (where the variables {q,p) are 
fixed and {q',p') are variables of integration with test functions a^q^p') G y{Q)): 

1 1 ^^^^ Akl p 

^^"^^2^"^^^^'^^'^''*'*'^" iV' ^ VdW - q - j^ - a + —{t - cTrev)]5{p' - p)} = 0. (40) 

fc=0 

The limit is performed as follows: ^ — )■ 0, a — ;■ 0, - — )■ 0, t = t{h), -{t — cTrev) -^ 2mD, 
hit — Y ) — ^ 0, where c G M, -D G [0, oo] and the numbers N' and a depend on c and N. If 
c is rational (and so, can be written as a reduced fraction c = ^), then N' = N for odd N 
and N' = Y for even N . Further, a = jj for N = 2 (mod 4) and a = otherwise. If c is 
irrational, then N' = 1, a = 0. The parameter a occurs in the definition of Ugp (see (J\) and 

a;, T_ = ^. 

The convergence in (^^ is uniform with respect to {q,p) on every subset of Q disjoint 
from some neighbourhood of the closed interval {p = 0} G Q. Moreover, if c is integer 
or half-integer, then the convergence is uniform on every subset of Q disjoint from some 
neighbourhoods of the points (±/,0). 

Since the meaning of the assertions of Theorem |2] is analogous to that in Theorem [H we 
repeat it only briefly here. Some complications appear, in particular, because we must take 
into account that the particle reflects from the walls and, hence, the sign of its momentum 
changes. Here ■^\{uqp,Uq'p^^t)\'^ is the probability density for a quantum particle to be in 
the state coqp at time t under the condition that it was in the state Uqipi at time (up to 
the norms ||wgp|P and Hw^yp, which tend to 1 in the sense of our limit on the sets under 
consideration; see Proposition [5] below. 

In the semiclassical limit ^, a, - — )■ a quantum particle in the state ojqp has a well-deflned 
position (equal to q) and momentum (equal to p),just like a classical particle. Therefore, 
we can say that, in the semiclassical limit, ^\{uqp,Uq'p'^t)\'^ is the probability density for a 
quantum particle in a box to be at the phase point {q,p) at time t under the condition that 
it was at the phase point {q',p') at time 0. 

As above, the case of c = and D = {— ^ 0) corresponds to the classical time scale 
Tcf. time is either flxed or increases slower than the rate of decrease of the collapse velocity 
(proportional to -).of the packet. Then the limit formula (HOj) takes the form 

1 T) 

lim[-—\{ujqp,Uq>p>^t)\^ - 6{q' -q + —t,p'- p)] = 0. 

Here S{q' — q + ^,p' — p) is the probability density for a classical particle in a box to be at 
the phase point (g,p) at time t under the condition that it was at the phase point {q',p') at 
time 0. Thus, in the semiclassical limit at time scale Td, we have classical dynamics: the 
quantum probability density of transition to {q,p) for a particle that was at {q',p') at time 
is equal to the corresponding classical probability density. 

The case c = 0, Z^ G (0, oo) (^ — )■ 2mD G (0, oo)) corresponds to the second time scale 
Tcoii- In this case we observe some spatial spread of the probability distribution. In contrast 
to the circle case, there is flattening not only of the spatial probability density but also of 
the probability of the signs of the momentum. 
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The case c = 0, D = oo (^ -^ oo) corresponds to a complete flattening of the spatial 
probability density (and the probability of the signs of the momentum): formula (140!) takes 

the form 

1 1 

lim ^ I Kp,Wg'p',t)r = -^^[S{p' -P) + S{p +p)]. (41) 

Therefore, we also associate this case with the time scale Tcou (corresponding to destruction 
of the localized wave packet). This yields a mathematical justification of the (asymptotic) 
flattening of the spatial probability density for a quantum particle in a box, which was 
previously known only from numerical experiments [El [7j. 

The case c ^ corresponds to the third time scale T^ev If c is irrational, then, as in the 
previous case, formula fllU]) reduces to fllTl) . that is, one observes a complete flattening of 
the spatial probability density. The case of rational c corresponds to a revival of the wave 
packet (full if A^' = 1, that is, c is integer or half-integer, and fractional otherwise). The case 
D > corresponds to the spread revived wave packets. In case D = oo, we again obtain a 
complete flattening of the spatial density of the distribution (141 p . 

The difference from the circle case regarding the condition p y^ for N ^ 1,2 (see the 
statement of Theorem |2]) is explained by a specific structure of fractional revivals of the states 
cuqo coming from the symmetry of these states (namely, from the property uj2i-qfi = — Wgo)- 

As in the circle case, since every irrational number can be approximated by rationals, we 
can (heuristically) say that the case of irrational c in Theorem [2] is a limiting case of rational 
c as A^ — )■ cxo: the distance between neighbouring terms in the sum of delta-functions in f HU]) 
tends to zero, and the sum of the delta-functions tends to the uniform distribution (in the 
weak sense). In other words, the cases of an irrational c and a very close rational c' are 
almost indistinguishable. 

Note that the distribution of the modulus of the momentum is preserved in all these 
limiting cases. 

Thus, we have traced the whole dynamics of a quantum wave packet in the infinite 
potential well. A well-localized wave packet eventually collapses with a complete flattening 
of the spatial density and the momentum sign. At some moments we see that copies of the 
initial packet simultaneously arise at several points of the well and then again eventually 
collapse with a complete flattening of the density. Since there are "more" irrationals than 
rationals, we can say that the particle most often stays in states whose spatial density of 
distribution is close to the uniform density. Thus, Theorem [2] completely describes the free 
quantum dynamics of a particle in the infinite potential well at all time scales the semiclassical 
limit. All stages are parametrized by the two real parameters c and D. 

Let us formulate a simplified version of Theorem [21 which deals only with principal time 
scales (classical motion, complete flattening, and exact revivals) without intermediate ones. 

Corollary 2. We have the following limit formulae in S^'{Q): 

1 77 

lim[^|(wgp, Wg'p',t)P - 5{q' -q + —t,p'-p)] = 
as h,a,- —^ 0, t = const; 
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2) 



1 11 



as h, a,- ^ 0, t ^ oo, — -)■ cx), 
we/ 
5; 



as well as h,a,^ ^ 0, t = cTrev — ^ 00, where c is irrational, T^ev = ~^: 



N'-l 



fc=0 

as /i, a, ^ — )■ 0, t = cTrev -^ 00, where c = ^ is rational and the numbers N' and a depend 
on N. 

The limits in cases 1), 2) and 3) for N = 1 or N = 2 (that is, for integer or half-integer 
c) are uniform with respect to {q,p) on every subset ofQ disjoint from some neighbourhoods 
of the points (±/,0). The convergence in case 3) for N ^ 1, N ^ 2 is uniform on every 
subset of Q, disjoint from some neighbourhood of the closed interval {p = 0} C Q. 

Proof of Theorem [H Since 

/2 

and [vqi-i^pi — vi^qi-pi) is an odd function, we obtain from the properties of G that 

{ujqp,UJq'p>^t)l = {Vq-l^p,Vq'-l^pi^t — Vl-qi -pi)2l 

(we will omit the subindex 21 in the remain part of the proof). We have 
1 



n 2vr/i 



[vq^i^p, Vqt-i^pt^t - vi-gi^-p>^t) I a{q,p)dq dp 

1 



^ t,,yVq~-i,p,Vq'p',t - v^q'-p'^t)\'^T ^ [a]{q' , p') dq' dp' . 
Furthermore, 

^1 (^^5-«,p, Vq'p'^t - l^-y'.-p'.t) P = ^ I {Vq~l,p, Vq'p',t)? + ^ | (^g-/,p, ^^-g',-p',t) ^ 

-Re{Vq_l^p,Vqi^p',t){V-q>-p>^t,Vq-l,p)- (42) 

Clearly, the last summand tends to zero uniformly on every subset of fi, disjoint from some 
neighbourhood of the closed interval {p = 0} (that is, one can find e such that |p| > e for 
all elements of the subset). 

Consider a weaker restriction: {q,p) lies in a subset of Q, disjoint from some neighbour- 
hoods of the points (±/,0). Since we already know that the summand tends to zero for 
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p 7^ 0, we consider the case p = 0. Let us expand the function T ^[a] into a Fourier series 
with respect to q' and represent it by a Fourier integral with respect to p': 






j=-oo 



aj(z/)e a'-' 



dq'dp 



Using formula fl^T]) and the method of the proof of Theorem [T|, we arrive at the formula 
1 



27rh 



{Vq-lfi, Vq'p>^t){V-q'-p',t, Vq-lfi)e'2l^ 



X^ exp 



4 
1 



2/" a2(4 + y) 



n 2 



4a2 



q — I — 2rl — - + uh 

4ml 



{q-l- 2rlf 



This expression tends to zero for g 7^ ±/ in the limiting cases corresponding to the case of 
rational c with A^ = 1 or A^ = 2 and, under this condition, the convergence is uniform on 
every interval [—1 + e, / — e], e > 0. That is. 



lim 



2vr/i././o 



>g-/,p, Vq\p'.t) {v-q',^p',t, Vq-i,p)T ^ [ct] {q , p) dq dp = 



(43) 



uniformly on {q,p) on every subset of fl disjoint from some neighbourhoods of the points 
(±/,0). 

Using equation fH2]) with zero last term. Theorem [H the definition of T~^ and the eveness 
property of T~^[a]{q',p') = T^^[a]{—q' , —p') of the image of T~^, we obtain for rational c 
that 

= lim / / ^Kwgp, w,y,t)|V(g',p') dq'dp' 

'J^l^ VD[q -q + l-j;^-a+-{t~—Trev)]T [a]{q,p)dq 

fc=0 "^"2' 

= lim // —^\{uqp,ujgy^t)\^(r{q',p') dq'dp' 



^ 2^ / VD[q - ^ - ^ - ^ + -(^ - ^^re.)]a(g ,p) rfg 

fc=0 "^"^ 
1 ^ /"' r / , 4fc/ p , M^ M / / X , / 

V -^^ j_i N' m IS 



where N' and a are defined as in the statement of the theorem. This proves the theorem in 
the case of rational c. 
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For an irrational c we have 
lim // ^;^\iujqp,(^g'p',t)\^criq',p')dqdp' 

'•21 1 r2l 



1 r^^ 1 z"^' 



4/ 

The uniform convergence in the two last formulae follows from that in flT9|) and the uniform 
convergence to zero of the expression ( l43l) or the last summand in f H2l) (depending on the 
limiting case under consideration). 

The theorem is proved. D 

Now let us prove an analogue of Proposition [2l 

Proposition 5. The norm ofcOgp tends to unity as h -^ 0, a ^ 0, - — )■ 0. The convergence 
is uniform on every subset of Q disjoint from some neighbourhoods of the points (±/,0). 

Proof. Using the properties of Ugp and 6, we have the following chain of equalities (see the 
beginning of the proof of Theorem [2]) : 

II l|2 _ / A _ /' ^ _ II l|2 / N 

W'^qpW — y^qp^^qpjl — [^q-l,p, ^q-l,p ~ ^l-q,-p)2l — \\^q-l,p\\ ~ \^q-l,p^'^l-q-pj2l- 

The first summand tends to 1 by Proposition [21 The second summand tends to zero for 
(g,p) G f2, {q,p) ^ (±/,0). Both limits are uniform on the sets indicated in the statement 
of the proposition. The proposition is proved. D 

We make three remarks on Theorem [2l 

Remark 6. The remarks at the end of Section [21 which were made in the circle case, remain 
valid here except for some inessential differences in the semiclassical dynamics of the mean 
position and mean momentum. First, the mean momentum is not preserved in the semi- 
classical limit in a box but periodically changes the sign (because of reflections in the walls) 
while its absolute value is preserved. Second, since the position of a particle in a box is well 
defined, we can use a familiar formula for the mean position: 

qt = x\uqp{x)f dx. 

Moreover, the centre of the packet is easily seen to move along the classical trajectory at the 
time scale Td. 

Remark 7. We have proved theorems on the semiclassical limit of quantum dynamics on 
a circle and on an interval (in a box). Of course, these results can easily be extended to 
the case of multidimensional domains that are Cartesian products of any number of circles 
and intervals: rectangles, rectangular parallelepipeds, tori, cylinders, and so on. It would 
be interesting to extend our results to the case of more general domains, including arbitrary 
domains in three-dimensional space, arbitrary compact manifolds, as well as to the case of 
interacting particles. 
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Remark 8. It is worthwhile to mention here the interesting results of M. V. Berry about the 
fractal images of the graph of \ip{x, t)p as a function of x and t. Here ip{x, t) is a solution of 
the Schrodinger equation for the infinite well (139|) with the uniform initial state ipi^x, 0) = -4=. 
In fact, the result is valid for a box of arbitrary shape and finite surface area in any finite 
dimensional space [SB]. 

Also, surprisingly, quantum revivals have an analogue in classical optics (Talbot interfer- 
ence) [371 [38] • 

4 The semiclassical limit of the dynamics for Husimi 
functions 

4.1 The case of free quantum dynamics on a circle 

A new (so-called functional) formulation of classical mechanics (or functional mechanics) 
was suggested in [TTl [H] (see also [191 1201 IHl [221 |23l [2l]). The basic concept of functional 
mechanics is not a material point or an individual trajectory, but the probability density 
function in a phase space. Accordingly, the fundamental dynamical equations are not the 
Newton (or, equivalently, Hamilton) equations but the Liouville equation (even in the case 
of one particle, not an ensemble). The Newton (Hamilton) equations become approximate 
equations for the mean values of the positions and momenta. Corrections to solutions of the 
Newton equations have been calculated in some particular cases [171 [T8l[2n [22] . 

The aim of passing to the new formulation of classical mechanics is to achieve the com- 
patibility of the reversible microscopic dynamics and the irreversible macroscopic dynamics. 
This problem is known as the irreversibility problem (or reversibility paradox). It is one of 
the most fundamental problems of mathematical physics. 

The dynamics of a material point is known to be reversible and recurrent while the 
dynamics of the density function satisfies the Liouville equation and has the so-called de- 
localization (collapse) property, which corresponds to irreversible behaviour. Therefore, if 
we adopt the description in terms of density functions not only for many-particle systems 
but even for one particle, then there is no contradiction between the micro- and macroscopic 
dynamics: both are irreversible (in some sense). 

A procedure of constructing the density function of a physical system from directly 
observable quantities (results of measurements) is described in [23]. The interaction of 
a system and the measuring instrument is studied from the point of view of functional 
mechanics in [23]. 

Now we want to know whether functional mechanics is preferable to Newtonian one from 
the quantum-mechanical point of view. To do this, we consider the semiclassical limits for 
quantum dynamics on a circle and in a box. In this subsection we take the case of the circle. 

There is a correspondence between quantum states (density operators) and classical states 
(distribution density functions), which was considered by Husimi in the case of coherent 
states on the whole axis [391 [30l [H] . Here we define an analogous correspondence for coherent 
states on a circle. 

Let p be a density operator (in other words, a quantum state) in L2{—l,l), that is, a 
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positive operator with unit trace. We associate with p the following function on the phase 
space Vt: 

piQ,P) = ^'^^^Mp = ^^{vgp,pvgp), (44) 

where Tr denotes the trace of an operator. If the density operator is a projector (p = P[ip], 
i^ G L2(— /,/)), which corresponds to the case of a pure quantum state, then formula 
takes the form 



Clearly, p{q,p) > and, by dl]), we have 

// p{(l^P)dqdp 



whence p{q,p) is a probability density function on the phase space (that is, a classical state). 
The correspondence (jSj) taking each quantum density operator to a classical density of prob- 
ability distribution is called the Husimi transform. It can also be expressed as smoothing the 
Wigner function of the quantum state with a Gaussian function [39lH0lBT] . The probability 
density function p is called the Husimi function of the operator p. 

Remark 9. The Husimi function is not the only way of mapping quantum density operators 
to classical probability density functions. One problem of the Husimi function is that its 
marginal distributions of position and of momentum does not coincide with the correspond- 
ing quantum- mechanical distributions. For example, if (for simplicity) p = P[ip] for some 
-0 G L2{—l,l), then the equality 

+ 00 

p{x,p)dp=\-^{x)\'^, (45) 

■oo 

in general, does not hold. Another way of mapping quantum density operators to classical 
probability density functions is the so called tomography map |l2l HS] HI] . A relation of the 
tomography map to the Husimi function is discovered in |15] . 

However, we take the Husimi transform for the following reason. Consider the positive 
operator-valued measure M given by the formula 



M{B) = ^jJ P[vgp]dqdp, 



(46) 



where B G flis a. Borel set. It can be regarded as an approximate simultaneous measurement 
of the position and momentum of a quantum particle on a circle (such measurements were 
introduced by J. von Neumann |16]). 

We would like to analyse the correspondence between classical and quantum dynamics on 
a circle. In classical mechanics, a simultaneous measurement of the position and momentum 
is allowed. Quantum-mechanically, this corresponds to an approximate measurement of the 
position and momentum like (1461) . 
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If a quantum particle is in the state p and we perform measurement (1461) . then the 
probabihty density function of the result is exactly p{q,p) given by dH]). 

Moreover, we will consider the Husimi function only in the semiclassical limit. In this 
limit, its marginal distributions coincide with the quantum- mechanical distributions (in par- 
ticular, dlS]) holds). 

More generally, a classical state a is a distributioncl. Suppose that a e ^'{K) (see Subsec- 
tion [3]T] for the definition of S>'{K)). The non- negativity of a is understood as (a, /) > for 
every non-negative test function /. The unit normalization of a distribution means that there 
is a limit lim (a, \p) = 1, where Xp, P > 0, is a family of functions such that X{x,p) = 1 

P— 5>00 

for \p\ < P, < X{x,p) < 1 for P < IpI < P + 1, and X{x,p) = for \p\ > P + 1. li a{q,p) is 
an ordinary function (a regular distribution), this condition expresses its integrability on Q 
and JJ^ a dqdp = 1. We denote the set of all non- negative and unit-normalized distributions 
from ^'{K) by ^[(K). 

We would like to know whether one can obtain any classical state as a semiclassical 
limit of the Husimi functions of some family of quantum states. This question can be posed 
rigorously as follows. Given a (distribution) probability density function a G ^[{K), can 
one find a family p^'^\ h > 0, oi density operators in L2(— /, /) such that p^^^ — )■ a in 2i'{K) as 
^ — )■ (where p'-^*-' is the Husimi function of the operator p*^^-*)? We note that the definition of 
the functions Vqp occurring in the Husimi transform involves a parameter a. As in Theorem [H 
we shall assume that a = a{h), a — )■ 0, - — ;■ (we have already noted that in this limiting 
case a quantum particle in the state Vgp has well-defined position and momentum). 

A further question concerns the dynamics. Let the quantum system (the density operator) 
evolve in time on a circle: pi = U^p^'^'^U^ , where t/f is defined by ([8]). We denote the 
Husimi function of Pt hy pi {q,p)- The evolution of the classical system (the probability 
density function) is given by the formula crt{(l,p) = <^{q — —t,p)- We may ask whether the 
correspondence between the family of quantum states p^'^\ h > 0, and the classical state 
a is preserved in time at various time scales (see Subsection 12.31) . Namely, is it true that 
PI — at ^ in ^'(K) as ^ — )■ 0, where t G M is an arbitrary fixed number (time scale Td)) 
or t — !■ oo (time scales Tcou and T^ev] see the various versions of the relation between the 
limits ^ — )■ and t — )■ oo in Theorem [T])? 

If the answers to these questions differ for a of the form <7{q,p) = 6{q — qo,P — Po) 
(which corresponds to an individual trajectory of a material point and, hence, to Newtonian 
mechanics) and for a G Li{Q) (which corresponds to a "bunch" of trajectories and, hence, 
to functional mechanics), then one can say that one formulation or the other is preferable. 

Remark 10. It is known that the asymptotic behaviour of at depends strongly on the initial 
density function a. If a{q,p) = 6{q — qo,p — po), then the motion is periodic and the period 
is — . If cr G Li(fi), then the spatial density asymptotically flattens with respect to q on 
a circle (in the sense of the so-called weak limit). Namely, the Kozlov's second theorem on 
diffusion (Theorem 2 in [13]) says that there is a limit 



lim a{q t,p)g{q,p) dqdp = — a{q',p)dq' 



giq,p)dqdp, (47) 



^In the most general case, a classical state is defined as a probability measure on the phase space [77] . 
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ii a,g G L2(fi). Formula (H7|) is easily seen to be valid for a G Li(fi), g G &{K). Thus, 

t— ^±00 



lim a(g - ^t,p) = 1 /" a(g',p) rfg' in ^'(iT), (48) 

:-5-±oo m 2f J_i 



if a G Li{Q). This is the sense in which the spatial density flattens. 
The following theorem answers the questions posed above. 

Theorem 3. 1) Let o G 2i'i{K). Then there is a family of density operators p^'^\ h > 0, in 
L2{—1, 1) such that the corresponding Husimi functions p^'^^ converge to a in ^'{K) as h —^ 0, 

' a 

2) Let, further, p\ = U^p^^'^Uf , and let pj he the Husimi function of pl . Then 



\im{'^!'\q,p,t) - a{q - —t,p)] = zn &'{K). (49) 



p_ 
m 

fi V n -A 4^( k\ fit 



Here the limit is performed as follows: /i — ?■ 0, a — ?■ 0, > 0, t = t{h), - — )■ 0. 

3) If a E Li{Q), then equality ( [TPP remains valid in the following limit: h ^ 0, a ^ 0, 
^^0,t^oo, ht^O, ^ ^ 2mD G (0, 00]. 

Proof, let a = <j{q,p) be a function from ^{K) (a particular case of distribution from 
S>'{K)) satisfying the conditions of non-negativity and unit normalization. Then we define 

P^''= [[cr{q',p')P[v,,,]p%. (50) 

The right-hand side depends on h because of the definition of the function Vg'pr. Moreover, 
we choose a in the definition of Vgp to be a function of h (that is, a = a{h)), such that a — )■ 
and - — )■ as /i — 7- 0. For example, we can take a = Cyfi. The integral is understood in 
the weak sense: p'-'*^ is an operator such that for all ^,x ^ L2{—1, 1) we have 

JjQ. W'^q'p'W JJq W^q'p'W 

We claim that p^'^^ is a density operator, that is, a positive operator with unit trace. Indeed, 
the positivity of p^^^ follows directly from the non-negativity of a. To prove that p'-^*-' is a 
trace class operator and to calculate its trace, let us take an arbitrary orthonormal basis Ui, 
i = 1, 2, . . ., in L2{—1, 1) and calculate 



f;(«.,p(%.) = j2ll ^i9',mu.,vgyr^^ 



ft °° do'dv' f f 

= (^{(l',P)y2\i^i^^g'p')\^Ti \H= // (r{q',p')dq'dp =1. (51) 

JJn fr^ Wvq'p'V JJn 

Here the third equation follows from the Parseval-Steklov identity. To justify interchang- 
ing the sum and the integral in the second equation, we note that since the sequence 
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Sj=i K""*)"^-?'?')!^' ^ ~ l)^,..., converges, it is bounded by some constant A. Then the 
integrand is majorized by the integrable function Aa and Lebesgue theorem enables us to 
pass to the hmit under the integral sign. Since the orthonormal basis was arbitrary, we have 
proved that p^'^^ is a density operator, as required. 
For t > 0, we have 






// (^{q',p')P[^Jq'p',t 



~{h)i X _ i / / I , ,.\i x|2 dq'dp' 



Pi iQ,P) = 7;-^ // (^iq\p)\iVqp,Vq'p',t)\' 



ZTChJjQ W'^q'p'W 

Using formula ( 14^ . Theorem [1] (the case c = 0, D = 0) and Proposition [21 we get 

lim[pl''^ (g, p) -a{q 1, p)] = 0, 

m 

uniformly on Q (and, hence, in ^'{K)), where the limit is performed as indicated in part 2) 
of the theorem. This proves parts 1), 2) of the theorem in the case when a G ^{K) (part 1) 
is obtained by putting t = 0). Parts 1), 2) for an arbitrary function a G ^\K) follow since 
y is dense in the set of locally integrable functions (that is, regular distributions), and the 
set of regular distributions is dense in the space 2> of all distributions [32] • 

Namely, to obtain an arbitrary distribution a G ^[{K) as a limiting case of a Husimi 
functions, we take a sequence of functions ar G ^{K), tendingto a in ^{K) as r — )■ oo and 
put 



P^^ = // ^rin)iq',p')P[v,v]T-^. (52) 

Here r{h) is a function tending to infinity as /i — )■ 0. Then the limit of the Husimi transforms 
is the desired distribution a. 

For a distribution of the form (j{q,p) = S{q — qo,p — Po), iqo,Po) £ ^ there is another 
regularization, which is more straightforward than the general one (just described). We put 

^(.) ^ P[^^l (53) 

1 1 ^qoPQ 1 1 

Then pf^ = p^^^ = P[^^^]. By Theorem [U (the case c = 0, D = 0) we have 
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pt\(l.P) = TT^, u^\iVgp,Vg^p^^t)\'^ ^ 6{q - qo t,p-po). (54) 

We now prove part 3) of the theorem. Suppose that o G Li{Vi) and p^'*-* is defined ( [50|) 
or (152|) . By Theorem [1] (the case c = 0, Z^ G (0, oo]) we have 

Y\Ta.[^'\q,p) - a{q 1 + q',p)^Diq') dq'] = 

J-i m 

uniformly on Q, where the limit is realized as indicated in part 3). On the other hand, we 
have the convergence (HHIl . Substituting this in the last formula, we see that 

\\m.pf^{q,p) = — j a{q,p)dq' (55) 
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in ^'{K). Comparing (155!) and (l48l) . we conclude that (149|) remains valid. The theorem is 
proved. D 

We note that formula (149|) does not hold if <j{q,p) = S{q — qo,p — Po), {qo,Po) & ^, and 
the limit is realized as in part 3) of Theorem [31 Indeed, define p^^^ by formula (153|1 . Then 






27Th\\Vq,p^^ 

lim[pj"^(g,p) - (pD{q - qo - —t)6{p - Po)] = 0, 
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where V9i5(g - go - ^i) 7^ <5(g - go - ^^) since D ^0. 

Formula fj49l) with a G -/^i(fi) also fails to hold if the limit is performed in the following 
way: ^ — )■ 0, a — )■ 0, ^ — )■ 0, t — )■ oo, t — ^T^ev — ^ 0, where ^ 7^ is a rational fraction and 
Trev is defined by formula ( IT^ . For example, if ^ = 1, then 

limpj''^(g,p) = a{q,p) 



while cr(g — ^t,p) satisfies formula 

Theorem [3] shows that an appropriate choice of a family of density operators enables us 
to obtain in the semiclassical limit any classical density function. Thus, in this sense, there 
are at least as many quantum states as classical states. This density function evolves in 
accordance with the laws of classical mechanics at time scale T^ independently of whether 
it is a delta-function or an integrable function. 

A difference appears at time scale Tcou- We have seen that if a is a delta- function, then the 
classical dynamics deviates from the quantum dynamics because a quantum packet collapses 
while a classical point-like particle remains the point-like particle for all time. Thus, the 
Newtonian mechanics holds only at time scale Td. From the other side, if a is an integrable 
function, then Theorem [3] shows that the classical dynamics remains valid from the quantum- 
mechanical point of view even at time scale Tcou since it exhibits a collapse of spatial density 
just as the quantum dynamics does. 

In both cases, the classical dynamics deviates from the quantum dynamics at time scale 
Trev because of the purely quantum effect of fractional revivals of packets. 

Thus, functional classical mechanics remains valid at a larger time scale than the New- 
tonian classical mechanics. Thus, it is preferable from the quantum-mechanical point of 
view. 

4.2 The case of quantum dynamics in a box 

Since the family Uqp, {q,p) G fi is a resolution of unity in L2(— /,/), the Husimi transform 
may be defined as 

p(g,p) = ^'^^^W/' = ^Kp.p^ap). 

where p is a density operator (a quantum state) in L2{—1, 1) and p{q,p) is a density function 
on the phase space Q (a classical state). If p = P[ip] for some tp G L2{—1, 1), then 

p{q,p) = ^'(^''^^p)!^- 
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The Husimi transform based on the functions Uqp is suitable for studying the correspondence 
between the classical and quantum mechanics for a particle in a box. We pose the same 
questions as in the case of a particle on a circle: is it true that every classical state can be 
obtained as a semiclassical limit of Husimi functions of quantum states, and at which time 
scale is this correspondence preserved for various classical states? 

A difference of the box case from the circle case is that u}±ifi = 0, whence we have 
p(±Z, 0) = 0. Therefore, we restrict ourselves to those density functions a which vanish in 
some neighbourhoods of the points (±/, 0). This restriction is physically inessential since the 
neighbourhoods can be arbitrarily small. 

Let the quantum system evolve in time in the box: pt = U^pUf , where U^ is the evolution 
operator in a box as defined by formula ( l38l) and p is an initial density operator. We again 
denote the Husimi function of pt by Pt{(l,p)- 

We can now prove an analogue of Theorem [3] for dynamics in a box. The non- negativity 
and unit normalization of distributions in ^'(f2) is defined as in !^'{K). We denote the set 
of all non- negative and unit-normalized distributions in ^'{Q) by !^[{Q). 

Theorem 4. 1) Let a G ^i(fi) and the support of a is disjoint from some neighbourhoods 
of the points (±/,0). Then there is a family of density operators p^'^\ h > 0, in L2{—l,l), 
such that their Husimi functions p^^^ converge to a in ^'{Vt) as /i — )■ 0, a — > 0, - — )■ 0; 

2) Let, further, pi = U^p^^^'U^ , pi and let pi he the Husimi function of p\ . Then 

\mi[^^\q,p,t) - a{q - —t,p)] = m ^'(fi). (56) 

m 

Here the limit is performed as follows: /?.—)■ 0, a —)■ 0, ^ — )■ 0, t = t{h), ^ -4 0; 

3) If a & Li{Q), then equality ([5^ remains valid in the following limit: /i — ;■ 0, a — )■ 0, 
^^0,t^oo, ht^O, ^ ^ 2mD e (0, oo]. 

The proof is similar to the proof of Theorem [3] but we now use Theorem [1] instead of 
Theorem [H In particular, if cr G S^{Q), then 



j j (^{(l\p)P[^q 



\U}q'p'\ 



and every distribution from S>'{VL) can be approximated by functions from S^{Q). Note that 
the integral on the right-hand side is well defined even though co±ifi = 0. This is because 
we assumed that cr vanishes in some neighbourhoods of the points (±/,0). But even if 
we take an arbitrary a G ^{fi), the points (±/, 0) are not poles for the integrand by the 
Cauchy-Bunyakovsky inequality: 



/ / n 1^ a{q,p)dqdp < / / a{q,p)dqdp. 



ZiTT ri/ J ./(> ll'-^g'T)' 



If a has the form a{q,p) = 6{q — qQ,p — po) for some {qo,Po) 7^ (i^jO), then we can 
propose a special approximating family of operators p*^^-' = P[t 



Thus, our conclusions for dynamics in a box are the same as for dynamics on a circle. 



,11 
fo 

The classical dynamics of a material point and the dynamics of a density function are both 

35 



obtained in the limit /i — )• at time scale Td. But the Newtonian fails to describe the collapse 
of the probability density and, hence, ceases to be valid at time scale Tcou- Functional 
mechanics, as well as quantum mechanics, exhibits the collapse of the probability density. 
But it does not describe the revival of wave packets, which is, thus, a purely quantum 
phenomenon having no analogue in classical mechanics. Hence, functional mechanics remains 
valid at time scale Tcou but not at time scale Trev We see that functional mechanics remains 
valid at a longer time scale than Newtonian one and, therefore, it is preferable from the 
quantum-mechanical point of view. 

Remark 11. In classical mechanics, analogues of the theorems on diffusion are proved for 
systems of interacting particles [H]. Therefore, it would be useful to generalize the obtained 
results to quantum dynamics of interacting particles on a circle and in a box, as mentioned 
in Remark [71 

5 Quantum dynamics in a box of size known with a 
random error 

We have so far assumed that the length 21 of the box is known with the infinite accuracy. But 
one of the basic postulates of functional mechanics says that the parameters of a physical 
system cannot be measured with the infinite accuracy. Hence, if we wish to predict the 
probability that the particle is in some place at a given time, we must replace the fixed 
value / G M by the density function /(/) of some probability distribution of the parameter 
I. In particular, one can construct this function from the results of measurements [23]. The 
dispersion of such a distribution may be very small, but it is always non-zero. Since T^ev 
depends on /, this breaks the exact periodicity of free quantum dynamics. Moreover, it turns 
out that the spatial density in a finite volume has a limit as t — ;■ oo. 

Theorem 5. Let ipi G L2{—l,l), I > 0, be any family of state vectors such that \iIji{x)\'^ is 
integrable with respect to (x,/) on [—1,1] x [0, oo), and let f{l) be a continuous probability 
distribution density supported on [0, oo) (that is, /(/) > 0, J^ f{l)dl = 1). Moreover, let 
f{l) = o{l^), e > 0, as I ^ Oj. Then, there is a limit 



hm P{x,t) = PUx)= r^^ 

t-s>±oo /n Zi 

where 



1 - ^iMy),My + 2X - 2/) +i;,{y-2x + 21)) 



f{l)dl, 



P{x,t)= / xi{x)mx,t)\^f{l)dx, 
Jo 

Xi{x) is the characteristic function of the interval [—1, 1], 4'i{x, t) = U^ipi{x), and the functions 
ipi are extended to the whole real line by the formula ipi{x + 2nl) = (— l)"-?/'i[(— l)"(x — 2nl)], 
n = ±l,±2,.... 

The function P{x,t) is the probability distribution density at time t taking account of 
the random error in the determination of /. The value of P{x,t) for a given x contains 
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contributions from the functions |-?/';(a;, t)p over all / such that x lies in the interval [—/,/]. 
This explains the factor xi{^) ^^ the expression for P{x,t). 

Let us analyse the limit probability distribution Poo{x). The term 



21 



f{l)dl 



corresponds to the uniform distribution: x;(a;)/2/ is the uniform distribution on the interval 
[— /, /], while the factor /(/) corresponds to the density of the probability that the half-length 
of the interval is /. The additional term 



Mx) 



Xi[x} 
4/ 



(My), My + 2a: - 21) +^,{y-2x + 21)) f {I) dl 



is, thus, a correction to the uniform distribution. Interestingly, it produces a dependence of 
the limiting final probability distribution on the initial one. We shall prove below that this 
correction tends to zero in the semiclassical limit. We also note that the limits as t — )■ +00 
and t — )■ —00, coincide, just as in the Kozlov's theorems. This reflects the time symmetry 
of quantum dynamics. 

Proof. Let us expand the functions tpi^x^t) into a Fourier series 



1 °° / 

i)i{x, t) = ^^ ak{l) sin ( 



nk . ,,\ f m'^hkH 



and substitute this in the expression for P{x,t): 
P{x,t)= I — ;— 2^ afcCnSm I — (x-/) I sm ( — (x-/) 



/ 



k,n=l 



X exp 



8mP 



f{l)dl. 



By the Riemann-Lebesgue lemma, all terms with k ^ n tend to zero as t — )> ±00. Hence, 

/ttA; 



hm P(x,t) = P^(x)= / M^^ 



t-5>± 



The theorem now follows from the formula 



|afcpsin2f^(x-/) ]f{l)dl. 



f:Krs,n^(^(.-o)4f: 

fe=l ^ ^ k=l 



\ak\ 

1 
2 



f Tik 
1-cosl — (x-/) 

1 - ]^{'iPi{y).My + 2a: - 21) + iji{y -2x + 21)) 



U 
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This result is rather unexpected since the quantum dynamics in a bounded domain is 
commonly regarded as being almost periodic [1] and, hence, has no limit for large values of 
time. 

We have already mentioned that the limit distribution differs from the uniform one. To 
estimate their difference A(x) in the semi classical limit, we first assume that ipi = ujqp is a 
coherent state (we omit the subscript /) and then represent an arbitrary state by an integral 
over coherent states. 

Proposition 6. Let ipi = cuqp for all I > 0. Consider the semiclassical limit h,a,- — > 0. 
Then A{x) tends to zero in the weak sense, that is, for every element a of the space ^(R) 
of rapidly decaying functions we have 

+00 

A{x)cx{x) dx — )■ 0. 

■oo 

Proof. Expanding a into the Fourier integral 

1 /"+00 



V27r y_oo 

and using the formulae ojqpiii + 2x — 21) = uj-q-2x-p{,y) and ujqp{y — 2x + 21) = U-q+2x~p{y), 
we get 

+ 00 -I /'+00 /»00 /' + 00 

A{x)a{x)dx = / dXa{X) dl f{l) {ujqp,uj-q-2x-p + (^-q+2x-p) dx. 

oo 4v27rf J^oo Jo J-oo 

Calculating the integral over x (the calculation is analogous to those in Section [3]), we see 
that this expression tends to zero. The proposition is proved. D 

Consider again arbitrary functions ipi G L2(— /,/), / > 0, and express them as integrals 
over coherent states (a corollary of fl36l) ): 



^ // 9iQ,P)^gpdqdp, 



Ztt 

where g{q,p) = {ujqp,ipi). We perform the semiclassical limit in the following way. As usual, 
let h., a, ^ tend to zero as parameters of the state Uqp., but let the function g remains constant. 
In other word, the dependence of ipi on h and a is such that the scalar product {uqp,ipi) 
coincides with the fixed function g{q,p). 

Corollary 3. Under this semiclassical limit of arbitrary functions ipi, we have A(s) — t- m 
the weak sense. 

The proof can be performed by a direct calculation. 

We recall that the weak convergence of probability density functions is sufficient from a 
physical point of view [131 HH [IS] • 

Thus, taking account of the inevitable random error in measuring the size of the box, 
we see that the probability distribution of the position of a particle in a box has a limit as 
t — > cx), and the limiting distribution is non-uniform and depends on the initial one. But in 
the semiclassical approximation, the limiting distribution becomes uniform. 
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Appendix 

Here we make a remark on the modular property of the theta-function 

+00 
9{z,t) = y. exp{— vrrA;^ + 27riA;z}, z, r G C, Rer>0, 
fc=— 00 
in a form applicable to the properties of quantum coherent states on an interval. 
We have the well-known modular relation 

e(-^\=^e^e{z,T). (57) 

\IT T J 

To see this, rewrite it in the form 

1 +00 f I \i-\ +°° 

1 v^ ii[z — nY\ sr-^ r 72 r^ -7 T 

—= 2_. 6xp < > = 2_^ expj— vrr/e + liiikz] 



(58) 



Express z G C in the form 2; = x + iy, x, y G M. The left-hand side of f lSS]) belongs to 
L2(0, 1) as a function of x for fixed y and r. Hence, it can be expanded into a Fourier series 
with respect to the orthonormal basis {e*'^'^'^, A; = 0,±l,±2,...}: 

4 E exp|- "^^ + ;^-"^ I = Y. ^^-'^'^- (59) 

* n=— 00 ^ ^ k=—oo 

We find the Fourier coefficients in the following way: 

, +00 „i r / , ■ \2 

1 Y^ / I 7t{x + ty — n^ 



* n=—oo ' ^ 

_ 1 /'+°° \ Tiix^iyf 

~ 7^ , 



a\. = -^= 2_, \ 6xp < mkx \ dx 

/'+°° f n(x + iyy ] 

/ exp < irckx > dx = exp{— tttA;^ + 27rik{iy)}. 



Substituting this expression in fl59|) . we obtain fl58|) . 

Thus, the modular property of the theta-function can be viewed as a rephrase of the 
Fourier series expansion of the function which is the sum of the Gaussian functions centred 
at the integer points. 
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